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I , Abstract 

p|J' We analyze radiation reaction for synchrotron radiation by com- 

K»" ■ puting, via a multipole expansion, the near field and derive from it the 

^ , Lorentz four-force, which we evaluate on the world-line of the charge. 

^ ' We find that the temporal component of the self four- force agrees with 

the radiated power, which one calculates in the radiation zone. This 
is the case for each mode in the multipole decomposition. We also 
find agreement with the Abraham-Lorentz-Dirac equation. 

1 Introduction and overview 

Radiation reaction is a century old problem of physics, which has been dealt 
with via many approaches during the years. (For reviews see, e.g., 0, ^ 
^.) Much of the understanding we currently have on radiation reaction 
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involves electromagnetic theory ||^, ^], but similar problems arise also in other 
branches of physics, e.g., in gravitation theory (and even in fields such as 
acoustics 1^). 

The usual approach for radiation reaction is based on arguments relating 
to the balance of quantities, which are constants of motion in the absence 
of radiation reaction, specifically the energy and angular momentum. Such 
balance arguments involve integration of the fiux of an otherwise conserved 
quantity over a boundary which normally consists of a distant sphere. This 
approach yields the energy and angular momentum radiated per orbit (i.e, 
the time-averaged energy and angular momentum fiuxes). Although such 
methods have been quite useful for the analysis of many problems involving 
radiation reaction, they are unsatisfactory for three main reasons. First, they 
apply only to periodic systems. Systems which involve radiation reaction are 
in a strict sense not periodic, because of the loss of energy through radia- 
tion. To illustrate the next point, consider a system comprising of a compact 
object in motion around a black hole, which is radiating gravitational radia- 
tion. Even if the system is approximated by a quasi-periodic model, assuming 
that the radiation reaction effects are small (specifically, that the time scale 
for inspiraling is much larger than the orbital period), eventually the orbit 
evolution due to radiation reaction will be fast. Recall now that the usual 
approach for radiation reaction yields only the time-averaged energy and 
angular momentum fiuxes. This leads us to the second problem with the 
usual approach: During the last stages of inspiral the evolution of the orbital 
motion is fast, such that the time-averaging methods would be inaccurate. 
The third problem occurs only for systems, which in the absence of radiation 
reaction involve non-additive constants of motion. The usual approach is 
based on a balance argument for conserved quantities. Specifically, knowl- 
edge of the radiated power at infinity implies that the system underwent a 
loss of energy at an equal rate, because of total energy conservation. How- 
ever, when non-additive constants of motion are involved, one can no longer 
add up the conserved quantity in the system and that part of the quantity 
which escaped to infinity, and then use the conservation law to deduce from 
the latter on the former. Such non-additive constants of motion occur for 
the very interesting case of motion in the spacetime of a spinning black hole, 
which is likely to be observed by future gravitational-wave observatories. 

An alternative approach is to calculate the radiation reaction directly, at 
the position of the system, rather than indirectly by balance arguments and 



calculation only of the far field in the radiation zone. The first such method 
was used by Dirac 0] for obtaining the Abraham-Lorentz-Dirac equation 
for an electric charge in arbitrary motion in Minkowski spacetime. In that 
approach one imposes local conservation laws on a tube surrounding the 
world-line of the particle, and integrates the conservation law across the 
tube, thus obtaining the equations of motion, including the radiation reaction 
effects. However, as is frequently the case with problems involving radiation 
reaction, the force which acts on the charge diverges. Dirac removed this 
divergence by using the so-called radiative Green's function, instead of the 
retarded Green's function. Specifically, one can write the retarded Green's 
function C*^* as the sum of two terms, namely 
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where G^'^'^ is the advanced Green's function. The first term is time symmet- 
ric. Namely, under the transformation t ^ —t the retarded field transforms 
to the advanced field and vice versa, such that the first term is transformed 
into itself. Because it is time symmetric, the first term does not include the 
radiative part of the field. Instead, it relates to the non-radiative Coulomb 
piece of the field, which is the source for the divergence. Dirac then showed 
that if one integrates only the second term, namely the radiative Green's 
function, one can obtain an expression for the field of the charge, which is 
finite on the world line of the charge. One can then use the usual Lorentz 
force formula to obtain the equation of motion which includes also the radi- 
ation reaction effects, i.e., the relativistic Abraham-Lorentz-Dirac equation. 
Dirac's prescription for regularization of the self force is quite successful for 
arbitrary motion in flat spacetime. However, because it involves the acausal 
advanced Green's function, it suffers from an inherent difficulty when one 
wishes to apply it for motion in a gravitational field: In curved spacetime 
the Green's function has support not only on the light cone, but also inside 
it. In particular, the instantaneous forces on the charge depend on the future 
history of the charge's motion, which is anti-causal [^]. Calculation of the self 
force which uses only the causal retarded field should therefore be preferred. 
Recently, Ori proposed a direct method for the calculation of self forces 
based on the retarded Green's function 0, ^. The method is based on the 
decomposition of the field and the self force into Fourier-harmonic modes. 
This approach has two important advantages: First, when one decomposes 



the field into modes one can separate the variables and obtain an ordinary 
differential equation rather than a partial one. Second, each mode of the 
self force turns out to be finite (even for point-like particles), such that di- 
vergences (when occur) are met only when one sums over all modes. The 
treatment of the individual modes is free from divergences. This approach has 
the advantages of being both local (such that it does not depend on balance 
arguments, which are useless for many interesting cases) and causal (such 
that it does not suffer from the acausalities of the radiative Green's func- 
tion). Of course, in general the sum over all modes of the self force diverges, 
and a regularization prescription is needed. Such a prescription has been re- 
cently proposed by Ori [T^. The availability of a regularization prescription 



strengthens the motivation for studying the self force using this approach. 
For many interesting cases the self force (or some of its components) is found 
to be finite without any need of a regularization procedure. These cases may 
serve as an exposition of the direct calculation of the self force, without the 
complications which are involved with regularization techniques. 

As an example for the application of this direct method for the calculation 
of radiation reaction we present in this paper an analysis of the radiation- 
reaction force for (electromagnetic) synchrotron radiation. Radiation reac- 
tion for this case turns out to be susceptible of an exact analytical calculation, 
i.e., one can sum analytically over all the Fourier-harmonic modes, and ob- 
tain a closed-form expression for the self force. In addition, the temporal 
component of the self force for synchrotron radiation, whose calculation is 
our main objective in this paper, turns out to be finite, without any need 
of regularization. In fact, also the angular components of the self force turn 
out to be finite: it is only the radial component which is divergent. (This 
divergence can be removed by simple procedures, and the renormalized radial 
self- force is zero.) The separation of the divergent component of the self force 



from the finite components is very useful |Tl|, as it allows the consideration 
of important aspects of the self force without getting involved in complicated 
regularization techniques. (In Section ^ we shall note in passing on the radial 
component of the self force.) 

Synchrotron radiation is an invaluable tool in the understanding of many 
physical phenomena, and is also very useful for many fields, such as astron- 
omy, chemistry, biology, medical science, and technology. For recent reviews 
see, e.g., [|12|, [l^, Q. The importance of synchrotron radiation stems from 



its unique properties, e.g., its spectral distribution in a broad range of wave- 



lengths, its relatively high intensity compared with other sources, the high 
degree of directionality of the radiation very close to the plane of the charge's 
orbit, and its well-defined polarization. 

The usual approach to radiation reaction from synchrotron radiation is 
to compute the far field in the radiation zone, then find the asymptotic 
expression for the Poynting vector and finally calculate the radiated power 
by integrating over a large sphere in the radiation zone. Then one expects the 
radiated power at infinity to equal the rate of loss of energy by the orbiting 
charge. Specifically, one assumes energy balance, and equates the radiated 
energy at infinity with the loss of energy by the charge. This approach 



appears in many standard textbooks, e.g., in Refs. [|1^, |16[. Here we develop 
a different approach, which is based on the near field rather than on the far 
field. Specifically, we consider the electromagnetic field which the orbiting 
charge produces, and find the interaction of this field with the very same 
charge that created it. Thereby, we construe the loss of energy by the charge 
as a self force or a back-reaction force. Our result here agrees with the well- 
known result, and also with the Abraham-Lorentz-Dirac equation [Q. Our 
method of calculation yields another important result, which is the multipole 
decomposition of the radiation. The radiation has been decomposed in the 
past into Fourier modes, but no previous multipole decomposition is known 
to us. 

We note that our entire discussion here is classical, and that we do not 
consider here quantum corrections. There are known cases where quantum 
corrections are extremely important, e.g., when there are extremely strong 
magnetic fields, such as those which exist in magnetars. The importance of 



quantum corrections is discussed in Ref. IT^- However, for most cases the 
quantum corrections for synchrotron radiation are unimportant, such that a 
purely classical discussion is useful. We stress that the agreement we find 
with the known expressions for synchrotron radiation or with the Abraham- 
Lorentz-Dirac equation is found just because we do not consider quantum 
correction. The Abraham-Lorentz-Dirac equation is purely classical, and 
does not include any quantum corrections. 

Another point which is worth emphasizing here is that our main result, 
the radiated power for synchrotron radiation, can be obtained very easily 



using conventional techniques, such as those which appear in Refs. [|1| 
Our derivation here is much more complicated from the technical point of 
view. However, from the conceptual point of view it is simple, as it is, in 



fact, only a direct calculation of the Lorentz force which uses only retarded 
fields. The technical complication is an unavoidable price which has to be 
paid in order to use this direct approach. 

We also note that the mathematical formalism which we use is more so- 
phisticated than the formalism which is usually used for the treatment of 
synchrotron radiation. The reason for this is that we consider the electro- 
magnetic field as derived from a four-vector potential. The decomposition 
of the four-vector potential into Fourier-harmonic modes is possible only in 
terms of vector spherical harmonics. Vector spherical harmonics are fre- 
quently used in standard texts (e.g., in Ref. []15[)- However, we preferred to 



use the conventions of indicial vector (and tensor) analysis, rather than the 



notation of []T3| . Because the analysis is done in curved coordinates, we shall 
frequently use concepts from Riemannian geometry (or pseudo-Riemannian, 
in a strict mathematical sense) . The reader is referred to any text on general 
relativity, e.g., Ref. IT 



The organization of this paper is as follows. We begin in Section ^ by 
decomposing the electromagnetic vector potential into vector spherical har- 
monics. In the Appendix A we briefly review the main properties of vector 
spherical harmonics which are important to us. More details can be found, 
e.g., in Ref. |T^. In Section ^ we compute the electromagnetic field from 



the vector potential, evaluate the field at the instantaneous position of the 
charge, and calculate the temporal component of the Lorentz four-force. We 
shall show that it is this component which is relevant for the radiated power. 
In Section ^ we use the more conventional approach and decompose the ra- 
diated power at infinity into spherical harmonics. This is done by evaluating 
the asymptotic far-field electromagnetic field, building the energy-momentum 
tensor from this field, and integrating the energy flux over a large sphere. 
This yields the radiated power at infinity directly, and we show that this 
far-field calculation agrees with the self force approach mode by mode. In 
Section ^ we collect the contributions from the various multipole moments, 
and expand the radiated power in a power series. This series can be summed 
over all orders, and the result we obtain for the total radiated power in Sec- 
tion]^ coincides with the well known result, which is obtained, e.g., in Refs. 
[p!5| , [TBI by completely different procedures. We demonstrate agreement with 
the Abraham-Lorentz-Dirac equation in Section |^, and discuss our method 
and results in Section Bl. 



2 The electromagnetic vector potential 

Let us consider a point-like electric charge q in circular orbit in Minkowski 
spacetime, such that the coordinates of the charge, in spherical coordinates, 
are r = ro, 9 = vr/2, = (pit), and t. The coordinates r,9,(f),t are inertial 
coordinates, e.g., coordinates of an inertial observer at r ^ tq. The proper 
time of the charge is r, such that t = 7r, 7 being the Lorentz factor. We 
assume the motion of the charge is uniform, 0(r) = 'yur. The angular 
velocity, as measured by a far away inertial observer, is u. We take the 
metric signature to be (— + + +). 



2.1 Temporal component of the vector potential 

We separate the four-vector potential A^ into a three-vector A and a scalar 
potential A*. Bold-faced quantities denote three vectors in any basis, namely 
either in a covariant, contravariant or in a normalized basis. (See Appendix 



A.) We shall use the Lorenz |2^ gauge throughout |^, and use units where 
the speed of light c = 1. The scalar potential satisfies the inhomogeneous 
wave equation 

DA* = -47rj*, (1) 

where the temporal component of the four-current density j* is given by 



j = pu = q'j dr 



r^ sin 6 



Here, p is the charge density, and u* is the temporal component of the four- 
velocity. We denote by □ the dAlembertian wave operator. Decomposition 
of j* into Fourier-harmonic modes yields 
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We denote ^''"(O) = e~'""'^y'™(cos^ = 0,0). In the Fourier decomposition 
we take the frequencies of the Fourier modes to be both positive and negative. 
We next decompose the scalar potential into modes according to 

/oo 1 

du-i,imAr)y''^{0,(t>)e-'^'. (3) 



Then, the left hand side (Ihs) of the wave equation becomes 
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A prime denotes differentiation with respect to the radial coordinate r. Sub- 
stituting into the wave equation (|l]) we find for the modes of the scalar 
potential: 
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We choose the solutions of the corresponding homogeneous equation to be 
such that one solution is regular at the origin (namely, at r = 0), and the 
other solution describes purely outgoing modes at infinity. These require- 
ments are equivalent to choosing a retarded solution, which is the physically 
relevant solution. These basic functions which solve the homogeneous equa- 
tion are 



and 
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ri 
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where j', h'- are the spherical Bessel and Hankel functions of the first kinds, 
respectively. Our choice of normalization (which is, of course, not determined 
from the homogeneous equations) is the same as in Ref . |]TP[ . The Wronskian 
determinant ||2^ of these two solutions is W = i sgn{u!) / (2,71) . We show in 
Appendix B that the solution for Eq. (^ can be written as a sum of two 
terms, each proportional to a Heaviside step function. Specifically, the Im 
mode of the scalar potential is given by 
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where 0(x) = 1 for a; > and Q{x) = for x < is the Heaviside step 
function, and where /(r) = —{ATTq/r)6{uj — mu;)P'™'(0). 



2.2 Spatial components of the vector potential 

The vector potential A satisfies, in the Lorenz gauge, an inhomogeneous 
wave equation 

DA = -47rj, (7) 

where j = pu is the three-current density of the sources for the field, and 
where u is the three- velocity. Decomposing the current density j into vector 
harmonics, we find 
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where Pj' '"(O) = e~*'"'^Yj ' ™(cos^ = 0, 0). We denote by a coordinate with 
an over-hat the coordinate in a normalized basis (see Appendix A), and by a 
star we denote complex conjugation. By Y- ' "^ we denote a component in a 
normalized basis of the pure-orbital vector spherical harmonics (see Appendix 
A). We next decompose the vector potential into modes according to 
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Then, the decomposition of the Ihs of the wave equation into modes yields 
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Each mode of the vector potential satisfies the equation 
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The basic solutions of the corresponding homogeneous equation are 
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such that the solution for the /', Im mode of the vector potential takes the 
form 
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where F{r) = —AnquSluj — mu) P^~'^"^ * (0) . 

We find it convenient to separate the electromagnetic field into electric 
and magnetic modes, which we later find (as indeed should be expected) to 
propagate on their own without mixing. The index /' of the pure-orbital 
vector harmonics (do not confuse the prime here with differentiation with 
respect to r) can take the values / — 1, /, or / + 1. After direct substitution of 
the explicit expressions for the various functions in Eq. (|12]) and summation 
over all /' modes, we find the Im mode of the magnetic part of the vector 
potential to be 

(13) 
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and, similarly, for the electric part 
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Here, 






llm 



I— 11m 



P'™(0) 



47rg?Ti u; ro 
.\/(/ + l)(2/ + l. 



AQ 



47rgr7i^ci}^ro ^ 
l{2l + 1 



a^ 



P'™(0) 



V^o 



-47rgmcjP'"(0) 



(20) 

(21) 

(22) 
(23) 



where dP^'^{0)/de = dY^"^ / de{e = 7r/2, = 0). The full vector potential 
is obtained by summing over the electric and the magnetic contributions, 
i.e., by A^ = ^A^^ + ^^ A^^ . By taking the complex conjugate of each 
component of the vector potential we summed, in fact, each m with its 
reciprocal negative counterpart. Specifically, because complex conjugation 
is equivalent to the transformation m -^ —m, adding the complex conjugate 
is equivalent to adding the — m counterpart to each term. Consequently, 
when we sum below over all modes, we need only to sum over positive values 
of m. 

We shall need below the following explicit expressions, which are easy to 



derive 23 : 
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2.3 The a; = modes 

We have not considered the possibihty that u = 0. This case corresponds 
to the magnetic number m = 0. In this case the terms proportional to 
u"^ on the Ihs of the wave equations (|^) and (|^) no longer appear. The 
independent solutions of the homogeneous equations are now ipQ = r^"*"^, 
ip^ = r~\ \E'q = r^"^^, and \E'^ = r~' . The Wronskian determinants are 
W = 21 + 1 and W = 21' + I, respectively. 

We find the vector potential for this case to be 

^ - 2/ + 1 ;^^ "^"^ "^"^ 2/ + 1 4^^ "^"^^ '^ ^^^^ 

2Z' + 1 H'+i ^ ° ^ 2/' + 1 r^+^ ^ °^ 

(27) 
It can be readily shown that this four- vector potential does not contribute to 
the temporal component of the self force or to the power carried off by the 
radiation, so we shall ignore it in the sequel. 



3 The self force 

We wish to present a calculation of the radiation-reaction force which is based 
on the near field of the charge. We first calculate the temporal component of 
the Lorentz force, namely, /* = qF^^Uji, where m^ is the four-velocity of the 
charge q, and F^p is the Maxwell field-strength tensor which is produced by 
the charge itself. This expression for /* is then to be evaluated on the world- 
line of the charge. The most convenient way to calculate the components 
of the Maxwell tensor is in covariant components, because then we have 
simply F^i, = A^^^ — A^^^, and we do not have to worry about covariant 
differentiation. (Although the metric is fiat, the spherical coordinates we 
use are curved, and therefore one would in general need to use covariant 
derivatives rather than partial derivatives.) Then, we shall change bases 
again to a normalized basis. We note that A^ = rsin^A?, Aq = rAg, A,. = 
Af, and At = Af. 

After calculating the temporal component of the self four-force we shall 
also calculate the radiated power, in order to demonstrate their compatibility. 

13 



We shall therefore be interested below in the radiated power, specifically in 
— -^ , where E is the energy of the charge. That is, -^ is the rate of 
change of the charge's energy, which is negative because the charge loses 
energy by radiating. The radiated power V is positive, and by conservation 
of energy V = -^. Specifically, -^ = JT'^ds^ = JT^'dsr, taking the 
surface of integration to be a large sphere, whose normal is radial. Then, 
dsr = T^ dQ, where dQ = sin 6 dd d(j). That is, we take for the radiated power 
at infinity — -^ = / T**" '™r^ dVt^ where T^y is the energy-momentum tensor 
of the electromagnetic field. However, the temporal component of the self 
four-force, by Newton's law, is the rate of change of the temporal component 
of the four-momentum, namely, the time derivative of the energy. In what 
follows we shall calculate both /* on the charge's world-line and the radiated 
power at infinity, and show that they agree. 



3.1 Magnetic modes 

The only non- vanishing components of the four velocity are the temporal and 
the azimuthal components. Namely, only Ut and u^ are not zero. Because the 
Maxwell tensor is skew-symmetric, it is clear that only F**^ can contribute to 
the temporal component of the Lorentz force. Calculation of this component 
yields 



^ t4> 
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When this is evaluated on the charge's world-line, we find that the field is 
continuous at tq. (We find that to be the case, in fact, for the entire temporal 
component of the self four- force.) We thus find that on the charge's world 
line 






mw 



charge 



l{l + l] 



-a 



Llm ■V' I — \ 



Qplm 



dd 



(0). 



(29) 
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Then, the contribution of the magnetic modes to the temporal component of 
the self force is given by 



M ri,lrn 



charge 



where u- 



charge 



q F u^ 



TqUJI- 



= 2g7^o r- — =ao 3 {^^^0}—- (0), 

charge Jl(l-U^^ OU 



l{l+l) 



(30) 



3.2 Electric modes 

Contributions to the electric modes come from /' = / + ! and /' = / — !. We 
denote the sum over these two values of /' by /. 



E JM 



£ .l+l,lm jt; .l-l,lm 



m 



(/ + 1)(2/ + 1) 



l+l,lm -1+1/ - \ 



ih ^^{mujr]e 



—imwt 



Y 



Im 



sin 9 



+ c.c 



m 



+ 



l{2l + 1 



m 



l—l,lTn -1—1/ — \ 



ih {■mujr)e 



-imuit 



Y 



Im 



(/ + 1)(2/ + 1) 



Cq J [muor) 



sin 9 



+ c.c 



Q{r -tq) 



sin 9 



+ c.c 



f^ l-l.lm -1-1/ - \ 

Oq J [muor) 



^1(21 + 1 



ih (mc<jro)e 



'imwt 



Y 



Im 



sin^ 



+ c.c 



0(ro - r) 



(31) 



Such that the t(f) covariant component of the Maxwell tensor is 



-.Im 



lira rp .Llm 



charge 






charge 



,2,-,, 



2a: 



l+l,lm m UJTq -l+l"^/^^-.^ \f=,lmt 



y(7TT)(27TTj 



j'+^ (mcjro)P'"(0) 



i2,-„ 



+ 2a'Q 



l-l,lm "^ ^^0 -/-iS/ _ N ^lmfn\ 



/(2/ + 1) 
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d^t 



+ 2^omf {mioro)P'"'{0) 

1 I'r i+l,iTn.i+l/ 



2nquj 



2" 



+ 



l — l,lm -l—l I - \ 



V9oj l^mcjroj 



(32) 



and the contribution of the electric modes to the temporal component of the 
self force is given by 



E f.*''™ 



charge 



^ tij),lm 

q F u^ 



charge 



2^A-A^,2f>lm2/n\ H' + -'-J 



(2^ + 1)2 



/ + 1 



j'+^(mcjro) 



1 _ _ n2 



(33) 



3.3 The total force 



Collecting the contributions of the magnetic and the electric modes we find 
that the Im mode of /* is 



f' 



tlm 



charge 



87r7(gm'cuVo)'/(/ + 1)T'™, 



where 



Here, 



r 



Im 



P''"(0)X 
2/ + 1 



/m,"| 2 



1 



dP' 



•Im 



^i^ = -/-i(mcjro 



ml{l + 1) de 
1 



-(O)j'(mc^ro) 



/ + ! 



j'+^(mc^ro). 



(34) 



(35) 



(36) 



In what follows we shall show that this expression coincides with the radiated 
power at infinity, and expand the force into a power series, which will enable 
us to find out how much radiation there is in each multipole moment. 

4 Mode decomposition of the radiated power 
at infinity 

We next expand the Maxwell tensor about r = cxo, and calculate the leading- 
order term in r~^ of the tr component of the energy-momentum tensor, which 
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is related to the radiated power. In each term of the energy-momentum tensor 
there is a product of two components of the Maxwell tensor. It can be shown 
that these products do not mix fields of different parities or components 
of different multipole moments. Namely, magnetic modes do not mix with 
electric modes (as should indeed be expected), and modes of different values 
of / or 772 do not mix. By a right arrow we denote that the quantity on the left 
hand side approaches the expression on the right hand side as r approaches 
infinity. We find 



^f[J , ^fJ -^ (_l)','+i4'""j'(mcuro)^=L=i^e^-^('--*) + c.c 



such that 



1 1 (9F'™ 

^l{l + 1) ^ -9^ 

(37) 



Fii F ^2ao' J^p^^-,3 {mur,)^^^^ + o.t. (38) 

We denote by o.t oscillatory terms, whose contribution to the total radiated 
power vanishes after integration over the sphere. Similarly, 

Fie , Fg- — >{-t)ao j {mujr)^===-——e ^ > + c.c (39) 

Jl{l + 1) " slllf^ 

such that 

Fie F ^2ao j^^^^-,J {n^r,)^^^ + o.t. (40) 

The total contribution of the magnetic modes to the energy-momentum ten- 
sor is thus 



Mrp f,lm _ J_ ( M r'™ M i7<^^.'™ I M r'™ M jj,6f,lm 

47r 



i A„ \ i<l> 'id I ■ \ ) 



We integrate this component of the energy-momentum tensor over a large 
sphere at infinity. Recalling that 

r I r)Y^"^ aylm * 2 _ \ 
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we find 



V 



M Im 



^^T*'''"^r2 dn 



1 l,lm2 ■I'^l - \ 



27r 



2_4 2 ^' 5P"(0)'.;2^ 



7(/ + l) a^ 
For the contribution of tlie electric modes we find 



^ t(f) 5 


^ <f>r 


and 






ET?lrn 
^ Or 


such that 





lAnqm^oo^ roP''"(0) ,„ F'"" 



2/ + 1 



-X' 



sin 6* 



^irnw{r—t) 



(42) 



+ C.C (43) 



—t 



,^,Anqm'u^r,P^^{Q)^i^dY 



Im 



2/ + 1 



X' 



dd 



,jimuj(r—t) 



+ C.C 






(44) 
(45) 



and we find 



V 



Elm 



^T'^'^'^r' dn 



,2_4,-4„2 H' + -'-J_pim2('Q\ 



57rg m cu Tq 



(2/ + 1) 



1 _ 1 _ _ "'2 

— — j'+^(mt^ro) - -j' ^(mu;ro) 



(46) 



The total radiated power for the Im modes is the sum of the magnetic and 
the electric contributions, which coincides with the expression above for the 
temporal component of the self four-force ( p4D up to a multiplicative factor 
of —7, which arises because the force is defined by /^ = dp^/ dr, whereas 
the radiated power is equated to dE/ dt . The minus sign results from the 
self force being the rate of change of the charge's energy which is negative, 
whereas the radiated power is positive. We stress that for the case of the 
radiated power we used the far-field expressions for the electromagnetic field, 
namely, the field very far from the charge, whereas for the case of the four 
force we calculated the near-field, namely, the interaction of the field which 
the charge itself creates with the very same charge. In both cases we used 
only retarded fields. 



1^ 



5 Power-series expansion of the radiation field 

In order to gain more insight into the expression we derived for the radiated 
power for each Im multipole moment, and to facihtate the summation over 
all modes, let us expand our expression in a power series. For each multipole 
mode we find a series expansion in cj^rg, such that only even powers of uro 
appear, as should be expected. This expansion is most easily obtained by 
using the standard expansions of the spherical Bessel functions for small 
values of their arguments. Next, we sum over all modes m for a given value 
of /. The power series expansion can be obtained for all values of /. For 
example, for the dipole electric mode (/ = 1) we find for the radiated power: 



^f = l<l'^'rl 



2_^ , 43 



i-i^'^« + is^'^« + «(' 



.-,6^6 



uj r. 



(47) 



This expression reflects the expansion of the radiated power in powers of the 
velocity squared. Our expressions for the radiated power in each Im mode 
of the field [Eqs. (^) and (^U])] are not restricted to slow motion, and are, 
in fact, fully relativistic. We find it convenient to expand in powers of the 
velocity (and thus obtain an expansion where the smallness parameter is the 
velocity) because this provides us with a very effective method to sum over 
all modes, which is translated into summation over all orders of the power 
expansion. When we sum over all orders in the expansion we are no longer 
restricted to the slow motion assumption, and we recover the fully relativistic 
nature of our analysis. 

Table |l] summarizes our results for the first three multipole moments for 
both electric and magnetic radiative modes, up to contribution of order uf'r^. 
(All higher modes contribute only at higher orders in the power-series expan- 
sion.) The first column of the table relates to the radiation mode, where Ex 
means the / multipole moment of electric radiation, and correspondingly Mi 
for magnetic modes. The second column describes the order of the contribu- 
tion. (All modes have a common factor of (2/3)g^ci;^rQ, which we suppress 
in this table.) The third column describes the contribution of the moment 
in question, and the fourth column says from what value of the magnetic 
number m the contribution comes. Recall that we, in fact, sum over the con- 
tributions of m and —m modes, such that the contribution of the m mode is 
just half of what Table |l| shows. 

We find that the dipole magnetic modes do not radiate, and that the 
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leading magnetic contribution comes from quadrupole radiation. This can 
be understood from the symmetry of the problem. Dipole magnetic radia- 
tion would correspond to a continuous current in a circular loop, say, with 
time-dependent current, which creates a magnetic moment, the second time 
derivative of which does not vanish. In our case the current results from 
the uniform motion of a discrete point-like charge, for which the second time 
derivative of the magnetic moment is zero. Therefore, the lowest moment for 
magnetic modes is quadrupole and not dipole as is the case for the electric 
modes. Because of the vanishing of the dipole magnetic modes, we find that 
the magnetic contribution is smaller than the electric contribution to the ra- 
diated power by a factor u'^Tq. Whereas the dipole magnetic contribution to 
the radiated power vanishes for all orders in cu^rg, higher multipole magnetic 
modes do not. In general, dipole modes should have contribution of order 
ui'^Tq, and for higher-order multipoles the leading contribution is smaller by 
a factor of uj'^^''~^h-Q . This is the case for both electric and magnetic 
modes, but because of symmetry the dipole magnetic mode does not con- 
tribute. Consequently, the leading magnetic contribution comes from the 
quadrupole modes. Another point which should be made is the contribution 
of the different m moments to a certain /. In the expressions we obtained for 
the radiated power, we see that electric modes have a multiplicative factor 

of P'™^(0), and magnetic modes have a factor of — ^^^ . Consequently, 
the contribution of electric modes vanishes whenever / — -m is odd, and the 
contribution of magnetic modes vanishes whenever l — m is even. We are now 
in a position to see why the magnetic modes are so weak: The dipole mode 
does not contribute because with / = 1, and bearing in mind that m = 
modes do not radiate, the only remaining m is m = 1, but then / — m = is 
even, and the contribution of the dipole magnetic mode vanishes. Therefore, 
the leading magnetic contribution comes from the quadrupole mode. How- 
ever, the m = 2 modes do not contribute to the quadrupole magnetic mode, 
because in that case / — m = is even, and it is just the m = 1 mode which 
can contribute. This results in the vanishing of the term of order ui^r^, and 
the leading magnetic contribution is only of order u^Tq. 
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6 The total radiated power 

We next sum all the contributions of the different mutlipole moments (the 
first few appear in Table ||). The summation is done in the following way. 
All the terms in Table |l| have a common factor of (2/3)g^c<)^rQ. The orders 
we shall refer to below are to powers of tu^rg beyong this common term. At 
order 1 only the electric dipole mode contributes, with coefficient 1. At order 
cu^rQ we have contributions from the electric dipole mode (with coefficient 
—2/5) and from the electric quadrupole mode (with coefficient 12/5). The 
sum of these contributions gives a total coefficient 2. At order u'^Tq we 
have contributions from the electric dipole mode (with coefficient 43/700), 
the electric quadrupole mode (with coefficient —16/7), the electric octupole 
mode (two terms, one with coefficient 729/140, and the other with coefficient 
1/2100), and the magnetic quadrupole mode (with coefficient 1/60). The 
sum of all these contributions gives a total coefficient of 3. The coefficients of 
higher-order terms can be found analogously. Collecting all the contributions 
up to order a)^ro, we obtain for the radiated power 

V = p'u'rl [l + 2u\l + ?>oj\i + O [u^t)] . (48) 

It can be shown (after a tedious calculation, or checking numerically) that 
for the nth order in the expansion inside the square brackets one finds the 
contribution to be nci}^"~^rQ"'~^, such that the general expression is 

9 oo 

V = -q^uj^rl £ ncD2"-V2"-2. (49) 
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This expression equals 



2 1 2 

T, ^ 2-4 2 -"- ^2 4-4 2 /cn\ 

^=^1"^ W, ^T^^ = o? 7 ^ '^o. (50) 

3 (1 — u^r-Q) 3 

which is the well-known expression for the radiated power of synchrotron 
radiation. 

We showed above that the Im mode of the radiated power corresponds to 
the Im mode of the temporal component of the self four-force. In fact, the 
only difference is an extra factor of —7, common to all multipole moments. 
Therefore, summing all the components of the self-force we find 

/ = -^g^Vrol (51) 
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However, the self force 

A dp^ dp^ dE 

as we indeed find. Here, p^ is the four-momentum. 

7 Comparison with the Abraham- Lorentz-Dirac 
equation 

We have obtained /* by calculating the Lorentz force acting on the charge, 
where the force originates from the electromagnetic field which the very same 
charge creates. On the other hand, the self force on the charge is given by 
the Abraham-Lorentz-Dirac (ALD) equation [E^, ffl, according to which 



where the operator D denotes covariant differentiation. (Recall that even 
though spacetime is flat, the coordinates are curved.) We shall next show 
that our result coincides with the ALD equation. Because the four-velocity 
is a constant four-vector, its partial derivative with respect to proper time 
vanishes. Therefore, the four-acceleration is given by 

"" ~ dr ~ dx-"" -^"/^"^ ' 

T^Q being the Christoffel symbols of the second kind. It can be readily 
verified that only the radial component of the four-acceleration does not 
vanish. Recall that only m* and u'^ are not zero. In spherical coordinates, 
the only non-zero component of the Christoffel symbols with lower indices 
t or is r^^ = — rsin^6'. (Note that F^^ = — sin6' cos 6', but the cosine 
vanishes on the charge's world-line, and therefore this component does not 
contribute to the four- acceleration.) Bearing in mind that sin^ = 1 on the 
world-line of the charge, we find that a*" = — ro7^ci}^, such that the squared 
four- acceleration is a^j^a^ = Vq'j'^uj'^. Substitution in the ALD equation (^) 
yields /ald = ~(2/3)g^7^ci}^rQ, namely, we find that 

/Ild = /*• (54) 
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Consequently, not only does our method for calculating the self force (by 
means of finding the Lorentz force acting on the charge by the field the very 
same charge produces) agree with the known expression for the radiated 
power at infinity, it also agrees with the Abraham-Lorentz-Dirac equation 

8 Discussion 

We presented a calculation for the temporal component of the self force for 
an electric charge in circular orbit in Minkowski spacetime. The calculation 
is based on decomposing the field which the charge creates into Fourier- 
harmonic modes, and finding the self force mode by mode. We derive the 
Lorentz force, and evaluate it on the world line of the charge. In the case 
of the temporal component it is possible to sum over all modes analytically. 
It turns out that this sum is finite, and is equal to the power carried by the 
emitted radiation to infinity. We also find agreement with the Abraham- 
Lorentz-Dirac equation. 

We note that even though it appears that the agreement between the self 
force we calculated and the radiated power at infinity is instantaneous, we 
in fact expect both quantities to agree in general only after averaging over 
a full cycle. We find in this case the agreement to be instantaneous because 
the problem is stationary. Consequently, it happens that the instantaneous 
value of the force is also the average over a cycle. In time- dependent problems 
we do not expect the agreement to be instantaneous. We still expect both 
quantities to agree over a full period, or after integration over all time if the 
motion is non-periodic. 

One might be surprised that the locally-derived self force turned out to 
be finite without any need of regularization. In general one would indeed 
expect the self force to diverge. However, it is reasonable to expect the self 
force to be the sum of a divergent piece and the finite, physical piece. The 
divergent piece is believed to be proportional to the acceleration. In our case, 
only the radial component of the acceleration is not zero. In particular, the 
temporal component of the acceleration vanishes, such that one can expect 
the divergent piece of the temporal component of the self force to vanish. 
This left us only with the finite piece, which is the physically important one. 

We briefly remark on the spatial components of the self four-force. From 
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symmetry, it is expected that the 9 component should vanish. This is in- 
deed what we find after a detailed calculation, similar to the calculation of 
the t component. The component is proportional to the t component. 
Specifically, 

f = qFKf, = qF^\ 

= -qF% = -—qF% 
euro ^ 

This component of the self force is the component which was calculated in 
Ref. |TI]. We note that one can ask the following question: What external 



force should be exerted on a charged particle, in order to keep it in uniform 
circular motion? In addition to the mechanical centripetal force, there is also 
another force, which balances the electromagnetic self force. The additional 
external force which should be applied is just the minus of the component 
of the self force we calculated here. (See also Ref. |TT|.) 



We also calculated the radial component of the self force using a similar 
approach. In that case, however, we were unable to sum over all modes 
analytically. This summation can still be done numerically. We find, that this 
sum diverges, although each mode is finite. The divergence can be removed, 
and the renormalized result agrees with the Abraham-Lorentz-Dirac result for 
the radial component, which is zero. We used two regularization methods for 
the radial component: First, we model the charge to be extended, and then 
we consider the limit of the size of the charge approaching zero. In a sense, 
this approach is in the spirit of the classical electron models of Abraham, 
Lorentz, Poincare, and others. The simplest model for a spatially-extended 
particle is a dumbbell model, where the axis of the dumbbell is radial. (We do 
not argue that this is a realistic model for a charge. One hopes, nevertheless, 
that the self force on a point-like charge would be independent of the model 
used in the limit of vanishing extension. In addition, the self force on a general 
extended object may be obtained by summing on all pairs of particles, where 
each pair is, in fact, a dumbbell.) After a simple mass renormalization (the 
bare mass is proportional to the acceleration of the dumbbell's center, and 
inversely proportional to the separation distance between the two edges of 
the dumbbell |2^), one indeed obtains the expected result of zero radial 
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force. A second renormalization method is based on the observation that the 
/ mode of the radial force on a point-hke charge approaches a constant as 
/ — i> oo. When one subtracts this asymptotic constant from each mode, the 
sum over all modes converges to the correct result. This latter regularization 



procedure is a direct application of Ori's prescription |T^, which is expected 



to be successful also for more complicated and interesting cases. 
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Appendix A: Vector spherical harmonics 

The electromagnetic vector potential is a four vector, and therefore it cannot 
be decomposed into the usual scalar spherical harmonics, but instead needs 
to be decomposed into vector harmonics. However, because the Minkowski 
spacetime is stationary (the metric functions are not functions of time), and 
thanks to the fact that the metric is diagonal, one can split the four vector 
potential into a scalar (the temporal component) and a three vector (the 
spatial components). The latter can then be decomposed into the usual 
vector spherical harmonics, whereas the former into the usual scalar spherical 
harmonics. Vector spherical harmonics are dealt with extensively in Ref. [|19| 



Here we describe just their properties which are the most relevant ones for 
our needs. We shall use the notation and convention of Ref. [IS]. We shall 



also give the relation of the vector spherical harmonic we use to those given 



by Ref. [|T|]. 



In order to decompose a general three-vector in flat spacetime, let us 
project it onto a tangent space to a sphere. Because this tangent space is 
two dimensional, we need two base vectors in order to span it. These base 
vectors are most conveniently chosen to be the gradient of the usual scalar 
harmonics, and its dual, which in two dimensions is also a vector [^]. These 
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are the Regge- Wheeler harmonics |^, defined by 



and 

where ai and hi are some normahzation factors, which we find exphcitly below. 
The tensor e^ can be represented by the matrix 

5 _ / -l/sin0 \ 
^'^ ~ l^ sin ^ y ' 

where the coordinates x" and x^ can take the values 6 or (p. The constants 
ai and hi can be found as follows. We take the Regge- Wheeler harmonics ipa 
and (pai and integrate their scalar product over the sphere. One finds 

^ ■ V-f *^f ^"^ dn = a*ail{l + l)5fi5„ 

and 



/0f*0iV^fi = &rw(/ + i)<5A 



and normalization yields ai = h = —r^ — . The metric q„h is the metric on 

the two-sphere, namely dT,"^ = Qabdx"' dx'^ = r^i^dO"^ + sin^6'(i0^). The cross 
product of ipa and 0a vanishes, namely 

i.lm * ilrh „ab 



Next, because the metric is diagonal, we can explicitly list the components 
of the Regge- Wheeler harmonics in a normalized basis [namely, in a basis 



defined by Vf,, = y |fi''^^|V^ (no summation over repeated indices) , where V^ 
is any four- vector written in a covariant basis]: 

1 (9F'™ 

"'* = 7I(7TT) ^^ 

/ 1 i^ vim 

^ y7(7+l)sin^ 
26 



«"^ ylm 



1 9F'™ 



We next denote ipa and 0a as F/' '™ and Fa^ '™, respectively, where x"' can 
take the values 9 or 0. In addition to the tangent space, we also have the 
radial direction, which is orthogonal to the tangent space. However, the pro- 
jection onto the radial direction would simply be just the scalar harmonics. 
Namely, we have y^^^™ = y'™. Y^, Y^, and Y^ are called collectively the 
pure-spin vector spherical harmonics. 

Finally, we find it most convenient for our needs to define the pure-orbital 
vector harmonics, because they are eigen-functions of the fiat spacetime 
Laplacian operator. These pure-orbital vector harmonics are built from the 
pure-spin harmonics in the following way: 

y/+l,/m _ (\/lY^^"^ — \/l + lY^^"^ 

yl-lM ^ ^ (y/TTlY^ '"" + VlY^ '" 

Namely, we find that y''''" is purely transverse, whereas the other two 
harmonics are two independent combinations of transverse and radial com- 
ponents. Table |^ gives the explicit expressions for the components of the 
pure-orbital vector spherical harmonics, in a normalized basis. 

For completeness, we also give the relation between the pure-spin vector 
spherical harmonics and the vector spherical harmonics which are given, e.g., 
in Ref. |]T3[ and involve the operator L = (l/?)r x V. First, we define 

vim T yim 



//(/ + 1) 

[This definition is identical to Eq. (16.43) of Ref. Il5|.] Then, 

-\/"S Im --vrlm 



27 



irxX' , 



and 



where f is a unit vector in the radial direction. 

It is useful to calculate the divergence of the pure-orbital vector harmon- 
ics. We find 



V • Y''-''" = -(/' - /)2 /' + hi' - / + 1) 



'/' + (/'-/)(/' -/ + 1)/21^,„ 



2/ + 1 r " 



where /' takes the values / — 1, /, or / + 1. The pure-spin harmonics satisfy 



2 
r 

B Im n 



VY 

Jl{l + 1) , 
r 
With these expressions for the divergence of the vector spherical harmonics, 
it can be readily shown that the Lorenz gauge condition is satisfied sepa- 
rately by each multipole moment of the field, by a direct substitution of 
the expressions we found for the components of the electromagnetic vector 
potential. 

Appendix B: Solution of a differential equation 
with a delta function source 

In this Appendix we shall show how one can find the solution to the inho- 
mogeneous second order ordinary differential equation 

g{x)y"{x) + q{x)y'{x) + h{x)y{x) = f{x)5{x - xq), 

when the two independent solutions u{x) and v{x) of the corresponding ho- 
mogeneous equation g{x)y"{x) + q{x)y\x) + h{x)y{x) = are known. By a 
prime we denote in this Appendix differentiation with respect to x. 
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We first assume that because the source has a delta function form, we 
shall have one functional form for the solution for x > Xq, and another form 
for X < Xq. Specifically, let us seek a solution of the form 

y{x) = Au{x)Q{x — xq) + Bv{x)Q{xq — x). 

Here, 0(x) is the Heaviside step function, which satisfies 0(x) = 1 for a; > 
and Q{x) = for a; < 0. The problem then is to find the constant coefficients 
A, B, such that we shall find the unique solution for the given inhomogeneous 
equation. Let us denote A{x) = AQ{x — Xq) and B{x) = BQ{xq — x). Then, 
the solution is of the form y{x) = A{x)u{x) + B{x)v{x). Differentiating 
A{x) and B{x) we find that A'{x) = A6{x — xq) and B'{x) = —B6{xo — x). 
Consequently, 

y'{x) = A{x)u{x) + A'{x)u{x) + B{x)v'{x) + B'{x)v{x). 

According to the variation of the parameters method, we now demand that 

A'{x)u{x) + B'{x)v{x) = 0. Substituting the expression for A'{x) and B'{x) 
we demand that A6{x—xo)u{x)—B6{xo—x)v{x) = 0, or [Au{xo)—Bv{xo)]S{x— 
Xq) = 0. Specifically, we find that B = Au{xo)/v{xo). 
Then, we find that 

y'{x) = A{x)u'{x) + B{x)v'{x) = Au\x)e{x - xq) + Bv'{x)Q{xo - x) 

and 

y"{x) = Au'{x)6{x—xq)+Au"{x)Q{x—Xo)~Bv'{x)6{xq—x)+Bv"{x)Q{xq—x). 

We next substitute these expressions in the differential equation. We 
use now the fact that the functions u{x) and v{x) solve the homogeneous 
equation to eliminate all the terms in the equation which are proportional to 
the Heaviside step function. We are thus left with 

g{x)Au'{x)6{x — xo) — g{x)Bv'{x)6{x — xq) — f{x)6{x — xq) = 0. 

The solution for this is 

Au'{xo) - Bv\xo) = fixo)/g{xo). 
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Recalling that B = Au{xo)/v{xo) we readily find that 

f{xo) v{xo) f{xo) v{xo) 



A 



g{xo) v{xo)u'{xo) - u{xq)v'{x(^) g{xo) W{xo) 



where W is the Wronskian of the two functions v and u. We thus find that 
^ = [fixo)/gixo)]v{xo)/W{xo), and B = [f{xo)/gixo)]u{xo)/W{xo), such 
that the solution for the differential equation is given by 

fjxo) v{xo) fjxo) ujxo) 

y{^) = ^-^TiTT^"^ ©^ ~ ^0 + ^-^T7T^^v{x)Q{xo - x). 
g{xo) W{xo) g{xo) W{xo) 
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Mode 


Order 


Contribution 


Value of m 


El 


1 


1 


1 




-2/5 


I 


,-,4^4 


43/700 


I 


Ml 


1 














,-,4^4 







E2 


1 







u^r^ 


12/5 


2 


uhi 


-16/7 


2 


M2 


1 







u^ri 







,-,4^4 


1/60 


1 


Es 


1 







u;X 







^'4 


729/140 


3 


1/2100 


1 


Ms 


1 














,-,4^4 








Table 1: The contribution to the radiated power of the first few multipole 
moments. First column: multipole moment in question. The capital letter is 
E, M according to whether the mode is electric or magnetic. The subscript 
is the corresponding /. Second column: the order of the contribution. Third 
column: the contribution at that order. Fourth column: the value of m 
wherefrom this contribution comes. Note that there is a factor of (2/3)g^u;^rQ 
common to all modes which is suppressed in this table, and that the E^ mode 
has contributions of the same order from two different values of m. 



33 



Component 


Expression 


yl + lMl 


V 21+1 ^ 


If 





yl-lM 

If. 


V 21+1 ^ 


yl+l.lm 




1 ey™ 


Va+i){2/+i) ^'^ 


•yl,lm 
8 


1 m i/'im 
^l{l+\) sine 


Yl-l,lm 

e 


1 9y"" 


V'(2«+l) «^ 


Yl+l,lm 

4> 


1 im yim 


v/a+l)(2/+l) ""^-^^ 


Yl,lm 


^ 1 9y'™ 


V'C+i) ^^ 


yl-ljm 

4> 


1 im i/'im 



Table 2: Pure-orbital vector harmonics in a normalized basis. The first 
column is the different components of the pure-orbital harmonics, and the 
second column is their value in terms of the scalar spherical harmonics. 
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